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Ph. We calculate the rate of production of + T and Z° + T in hadron 



colliders. We find the cross sections for + T and Z° + T to be roughly 



^ . 0.45 pb and 0.15 pb at the Tevatron and roughly 4 pb and 2 pb at the LHC. 

The dominant production mechanism involves the binding of a color-octet bb 
pair into a P-wave bottomonium state which subsequently decays into T. The 
purely leptonic decay modes of T, W^, and Z° provide signatures with small 
backgrounds. These events may be observable in Run II at the Tevatron, and 
they should certainly be observable at the LHC. 
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In high energy collider experiments, the particles with the cleanest signatures are those 
with purely leptonic decay modes. They include the charged weak gauge bosons W^, which 
decay via W — > £v where £ is an electron or muon, and the neutral weak gauge boson Z°, 
which decays via Z° — > £ + £~ . They also include the J = 1 quarkonium states, such as 
the charmonium state J/ip and the bottomonium state T, which decay into £ + £~ . Purely 
leptonic decay modes are particularly useful in hadron colliders, because they provide an 
enormous suppression of the background. Events involving particles with such decay modes 
can be employed both as a probe of the production process and as a lamppost under which 
to search for new physics. In this letter, we calculate the rate in the Standard Model for the 
associated production of the bottomonium state T with W or Z° at the Tevatron proton- 
antiproton collider and at the LHC proton-proton collider. We find that the dominant 
production mechanism for both W + T and Z° + T involves the binding of a color-octet 
bb pair into a P-wave bottomonium state which subsequently decays into T. The purely 
leptonic decays of the T, W^, and Z° provide clean signatures for these events. The cross 
sections for + T and Z° + T may be too small for these events to be observed at the 
Tevatron, but they should certainly be seen at the LHC. 

The associated production of W ± or Z° and a J/ip has been studied previously, but 
only for J/^'s with large transverse momentum, where the process is dominated by gluon 
fragmentation into J/ip The differential cross sections for W ± + J/ip and Z° + J/ip peak 
at small transverse momentum of the J/ip. Unfortunately, if the J/i/j is produced with small 
transverse momentum, one or more of the leptons from its decay is likely to have insufficient 
energy to be identified. Because the mass of the T is larger than that of the J/ip by a factor 
of 3, the cross sections for W + T and Z° + T are much smaller than those for W ± + J/ip 
and Z° + J/ip. However the larger mass of the T makes it possible to observe the leptons 
from its decay, even if the T is produced at small transverse momentum. It is therefore 
possible to measure the total cross sections for W ± + T and Z° + T. 

The production of T requires the creation of a bb pair at short distances of order l/m& or 
smaller. The bb pair subsequently forms a color-singlet bound state, which can either be T 
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or a higher bottomonium state that decays into T. The probability of binding depends on 
the bottomonium state and on the color and angular momentum quantum numbers of the bb 
pair. In the NRQCD factorization approach to inclusive quarkonium production 0, these 
probabilities are parameterized by nonperturbative matrix elements in nonrelativistic QCD 
(NRQCD). The matrix elements scale in a definite way with the typical relative velocity v 
of the bottom quark in bottomonium. Its value is roughly v 2 m 1/10. The NRQCD scaling 
rules can be used to estimate the order of magnitude of the binding probability, with every 
factor of v 2 corresponding to suppression by an order of magnitude. 

We first consider the direct production of T. The probability of forming an T is largest 
for the bbi( 3 Si) state, which is a color-singlet configuration with angular-momentum quan- 
tum numbers 3 Si. It can be parameterized by the NRQCD matrix element {Oj( 3 Si)) , which 
is proportional to the square of the wavefunction at the origin and can be extracted phe- 
nomenologically from the leptonic width of the T: (Oj( 3 Si)) = 7 GeV 3 . According to the 
velocity-scaling rules of NRQCD, the bb configurations with the next largest probabilities 
for binding into an T are the color-octet states bb 8 ( 3 Si), ^sO^o), and bb s ( 3 Pj) 0. Their 
binding probabilities are suppressed by v 4 relative to 5&i( 3 Si). If the cross section for color- 
octet bb pairs was comparable to that for bbi( 3 Si), then the color-octet contributions to 
the cross section for direct T production would be suppressed by two orders of magnitudes. 
However in hadron colliders, the cross section for color-octet bb pairs is much larger than 
that for bbi( 3 Si), and this can compensate for the smaller binding probability. The cross 
section for &6 8 ( 3 Si) is particularly large, because it can be produced by the decay of an 
off-shell gluon. There are as yet no reliable phenomenological determinations of the matrix 
elements (Oj( 3 Si)), (OJ^Sq)), and (Oj( 3 Po))- In their analysis of T production at the 
Tevatron ||, Cho and Leibovich used the value (Oj( 3 Si)) ~ 0.006 GeV 3 , which they ob- 
tained by taking a phenomenological value for (Og^ ( 3 S\)} with large error bars and scaling 
it by (mbVb) 3 /(m c v c ) 3 . They extracted a value for the linear combination (Oj (^Sq)) /5 + 
(Oj ( 3 ; P )) / ' m l by fitting CDF data for bottomonium production ||. They obtained a cen- 
tral value 0.008 GeV 3 with a statistical uncertainty that was greater than 100%. Thus the 
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uncertainties in the color-octet matrix elements are probably at least an order of magnitude. 

We next consider the indirect production of T by the binding of the bb pair into a higher 
bottomonium state that subsequently decays into T. The higher states include T(2S) and 
T(35), which decay into T with branching fractions of about 31% and 13%, respectively. 
Since these states have the same quantum numbers J PC = 1 as T, their NRQCD matrix 
elements scale in the same way with v and will therefore be comparable in magnitude to 
those for T. For example, the binding probabilities of &&i( 3 Si) into T(2S) and T(3S) are 
both about 40% smaller than for T. The effect of the feeddown from higher bottomonium 
states on the production of T from bbi( 3 Si) can be taken into account in the NRQCD 
factorization formulas by replacing (Oj^Si)) by 

J2B(H -> TX)(Of( 3 S 1 )) = 8GeV 3 , (1) 

H 

where we have added the contributions from T (with B = 1), T(2S), and T(3S f ). Including 
the feeddown from T(2S) and T(3S) therefore increases the cross section by about 20%. 

Bottomonium states with different J PC quantum numbers can have a more dramatic 
effect on the cross section, because they receive contributions from parton processes with 
different selection rules. The only process whose binding probability is of the same order 
in v as 6&i( 3 Si) into T(nS) is bbi^So) into rjb{nS). By the heavy-quark spin symmetry 
of NRQCD, the matrix element (Oj fc(,lSf) ( 3 S , i)> differs from that of (Oj {nS) (}So)) only by a 
spin factor of 1/3. Including the contributions from r] b (2S) and r] b (3S) and assuming that 
B(i]h(nS) — > TX) is comparable to the branching fraction for hadronic transitions from 
T(nS) to T, we obtain the estimate Eh B(H -> TX)(Of ( l S )) « 0.3 GeV 3 . 

The binding probabilities for P-wave bottomonium states are suppressed by v 2 relative 
to T. For Xbj{nP), the bb configurations with binding probabilities suppressed by only v 2 
are bbi( 3 Pj) and 66g( 3 5'i). By heavy-quark spin symmetry, the color-singlet matrix elements 
(Of j(nP) ( 3 P/)) differ for J = 0,1,2 only by a spin factor of 2 J + 1. Their values can be 
determined from the wavefunction in nonrelativistic potential models. We adopt the values 
used in the analysis of Ref. ||. Weighting the matrix elements for Xbj(lP) and x&j(2P) 
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by their measured branching fractions into T, we obtain J2h B(H — > TX)(0^ ( 3 Pj)) /ml ~ 
(0.003,0.2,0.2) GeV 3 for J = (0, 1, 2). We have included a factor of 1/mg in the P-wave ma- 
trix elements so that they have the same dimensions as S-wave matrix elements. The relative 
magnitudes of the matrix elements can then be interpreted as relative binding probabilities. 
The color-octet matrix elements (Og bj( ' nP \ 3 Si)) differ for J = 0, 1,2 only by a spin factor 
of 2 J + 1. In the analysis of T production at the Tevatron by Cho and Leibovich ||, these 
matrix elements were determined by fitting the CDF data on bottomonium production [|J . 
Weighting the matrix elements by their measured branching fractions into T, we obtain 

Y,B(H ^TX)(Og( 3 S 1 )) « 0.4 GeV 3 . (2) 

H 

We have included the contributions from Xbj(lP) with J = 1,2 and from Xbj{2P) with 
J = 0, 1, 2. The IP states account for about 70% of the total. The value (0) is more than 
an order of magnitude larger than the estimate by Cho and Leibovich of the matrix element 
(Oj( 3 Si)) that contributes to direct T production. The uncertainty in (0) is probably 
at least a factor of 3, because the analysis of T production has substantial theoretical 
uncertainties. The weighted matrix element in (^) is smaller than the corresponding color- 
singlet term in (|I|) by an order of magnitude. This smaller binding probability is more than 
compensated by the much larger cross section for bbg^Si) in hadron colliders. Thus the 
binding of &6 8 ( 3 >S'i) into Xbj followed by its decay into T gives a large contribution to the 
total T cross section. 

We now consider the production of W ± + T in a hadron collider. This process proceeds 
through parton collisions that produce final states including W+bb. The relative importance 
of the various contributions depends on the magnitude of the NRQCD matrix elements and 
on the size of the parton cross sections. The only parton cross sections that can be calculated 
relatively easily are those that can be taken into account via 2 — ^ 3 parton processes of the 
form ij — > W + bb. Fortunately, they account for the largest contributions. The 2-^3 
parton processes that produce W + + bb are 

• Wl: ud, cs — > W + + bbs^Si) via the decay of a virtual gluon into bb, 
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• W2: cb — > W + + bbi^Si, 1 S , 3 Pj) involving the exchange of a virtual gluon, 

• W3: ud, cs — * W + + bbi( 3 S\) via the decay of a virtual photon into fofo. 

We have omitted processes that are suppressed by a tiny Kobayashi-M askawa factor. Process 
Wl is the only color-octet contribution with a cross section of order a 2 s a w . All other color- 
octet contributions are suppressed by an additional factor of a s . Color-singlet contributions 
involving QCD interactions are suppressed by an additional factor of a 2 s . One such process, 
gg -> W + + bb^Si) + foe, which has a cross section of order dominated by the 

splitting of the colliding gluons into a collinear bb pair and a collinear cc pair. It can be 
taken into account through process W2. We expect this process to provide a reasonable 
estimate of the color-singlet contributions involving QCD interactions, because it is the only 
such process that is enhanced by a factor of \og(M w /mb) \og(M w /m c ). Process W3 is a 
color-singlet contribution that involves only electroweak interactions and has a cross section 
of order a 2 a w . 

We proceed to calculate the cross section for production of W ± + T at the Tevatron. We 
use the leading order GRV-94 parton distributions, and set the factorization and renormal- 
ization scales equal. We neglect the masses of the u, d, s, and c quarks, and we use the value 

4.7 GeV for the b quark mass. The cross section in pp collisions at center-of-mass energy 

1.8 TeV is 

Y>(W± + T) = 1.6 ft XmWL + 450 fb SW^)>. (3) 
± 8 GeV 3 0.4 GeV 3 V ; 

We have omitted the contributions from the color-singlet 1 S'o and 3 Pj matrix elements, be- 
cause they are several orders of magnitude smaller than the color-singlet 3 Si term. The 
color-singlet 3 Si term is completely dominated by the purely electroweak process W3. How- 
ever it is overwhelmed by the color-octet 3 Si contribution from process Wl. While the 
color-singlet matrix element is accurately normalized by the leptonic decay width of the T, 
there is a large uncertainty in the normalization of the color-octet matrix element. The the- 
oretical error due to higher order perturbative corrections can be estimated by varying the 
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factorization and renormalization scales by a factor of 2 and is roughly 35%. The theoretical 
error due to the uncertainty in the b quark mass can be estimated by varying m& by 0.25 
GeV and is roughly 25%. 

To obtain the cross section in the purely leptonic decay channels, we multiply by the 
branching fractions 5% for T — > £ + £~ and 22% for W — > £v. Given the integrated luminosity 
of 110 pb -1 in Run I of the Tevatron, there should be about 0.5 events in the leptonic 
decay channels. To estimate the number of events that would actually be detected, we 
must multiply by the detector acceptances and efficiencies. The product of the acceptances 
and efficiencies for the CDF detector was about 26% for W — > ev || and about 12% for 
T — > [|]]. Thus after allowing for acceptances and efficiencies, the cross section is 

about two orders of magnitude too small to be observable in the data from Run I. 

We next consider the production of Z° + T in a hadron collider. The 2^3 parton 
processes that produce Z° + bb are 

• Zl: uu,dd,ss,cc^ Z° + bb 8 ( 3 S 1 , 1 S , 3 Pj), 

• Z2: gg^ZV + bhfS^So'Pj), 

• Z3: gg^Z^ + bb^S^S^Pj), 

• Z4: bb -»• Z° + bb 1 ( 3 S 1 , 1 So, 3 Pj), 

• Z5: uu, dd, ss, cc — ^ Z° + bbi^Si) via the decay of a virtual photon into bb. 

Processes Zl and Z2 are color-octet contributions with cross sections of order a 2 s a w . Process 
Z3 is the only color-singlet contribution with a cross section of order a 2 s a w . Another color- 
singlet contribution gg — > Z° + 66 1 ( 3 5'i) + bb, which has a cross section of order a^a w , can be 
taken into account through process Z4. Finally, process Z5 is an electroweak color-singlet 
contribution with a cross section of order a 2 a w . 

The cross section for production of Z° + T in pp collisions at the Tevatron with center- 
of-mass energy 1.8 TeV is 
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a(Z° + T) = 8.5 fb 



8 GeV 3 



+ 150 fb 



0.4 GeV 3 



(4) 



We have calculated all contributions to direct T production with binding probabilities of 
order v 4 or larger and all indirect contributions involving binding probabilities of order v 2 or 
larger. All those terms that are not shown explicitly in (f|) give contributions that are more 
than an order of magnitude smaller than the color-singlet 3 Si term. The color-singlet 3 Si 
contribution is dominated by processes Z3 and Z5, which contribute about 85% and 15%, 
respectively. However it is overwhelmed by the color-octet 3 Si term, which is dominated 
by process Zl. While the color-singlet matrix element is accurately normalized, there is a 
large uncertainty in the normalization of the color-octet matrix element. We estimate the 
theoretical errors from higher order perturbative corrections and from the uncertainty in the 
b quark mass to be roughly 35% and 25%, respectively. 

To obtain the cross section in the purely leptonic decay channels, we must multiply by 
the branching fractions 5% for T — > £ + £~ and 7% for Z° — > £ + £~. Given the integrated 
luminosity of 110 pb -1 in Run I of the Tevatron, there should be about 0.06 events in the 
leptonic decay channels. The product of the acceptances and efficiencies of the CDF detector 
was about 30% for Z° — > e + e~ || and about 12% for T — > ji + //~ ||. Thus after allowing 
for acceptances and efficiencies, the cross section is several orders of magnitude too small to 
be observable in the data from Run I. 

In Run II of the Tevatron, the increase of the center-of-mass energy to 2.0 TeV will 
increase the cross sections for W ± + T and Z° + T by about 10%. With an integrated 
luminosity of 2000 pb _1 , the number of events in the purely leptonic decay channels should 
be about 10 for W ± + T and about 2 for Z° + T. The upgrades of the CDF and DO detectors 
should increase the acceptances and efficiencies for observing these events. Given that the 
uncertainties in the cross sections are at least a factor of 3, there is a possibility that these 
events could be observed in Run II. An observation of these events at a much larger rate than 
predicted could be evidence for a heavy particle that has a substantial branching fraction 
into W ± + bb or Z° + 66. 
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The cross sections for production of W ± + T and Z° + T in pp collisions at the LHC 
with center-of-mass energy 14 TeV is 



J2^(W ± + T) = 10 fb 



8 GeV 3 



+ 4000 fb 



a(Z° + T) = 5OOfb EjB(01 ff l)) - 1300 11. 



0.4 GeV 3 ' 
J2B(0 8 ( 3 S 1 )) 



(5) 
(6) 



8 GeV 3 0.4 GeV 3 

With an integrated luminosity of 10 fb -1 , the number of events in the purely leptonic decay 
channels should be about 440 for + T and about 70 for Z° + T. Even after allowing for 
detector acceptances and efficiencies, there should be enough events to make these processes 
observable. In Figure [l], we plot the invariant mass distributions da/dMwr (summed over 
H /± ) and da/dMzr- They peak at only a few GeV above the thresholds of 89.8 GeV for 
W ± + T and 100.6 GeV for Z° + T. 
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FIG. 1. The invariant mass distributions do / dMy/T (upper curve) and da/dMzr (lower curve) 



for the production of + T and Z° + T in pp collisions at center-of-mass energy 14 TeV. 
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The cross sections for + T and Z° + T are both dominated by the (Og( 3 Si)) 
term. The uncertainty in the value of B(0$( 3 Si)) leads to a substantial uncertainty 
in the normalizations of the cross sections. However this uncertainty cancels in the ratio 
J2± &(W ± + T)/cr(Z° + T), which is predicted to be about 3 at the Tevatron and about 2 
at the LHC If a significant deviation from this prediction is observed at the LHC, it might 
be evidence for an additional contribution from a heavy particle that decays into + bb 
or Z° + bb. 

The production of W + T and Z° + T provide lampposts under which one can look 
for new physics. The most promising possibility is to search for a charged Higgs via the 
decay H + — > W + + T. The decay rate of the charged Higgs into W + bb is enhanced by the 
Yukawa coupling of the Higgs to a virtual top quark. If the mass of the charged Higgs is 
in the range between 140 GeV and tb threshold and if the Higgs mixing parameter tan (3 is 
small, then W + bb may be the largest single decay mode ||. The decay rate of the charged 
Higgs into W + T was first calculated by Grifols, Gunion, and Mendez J/J . For small tan /3, 
the branching fraction B(H + — > W + + T) ranges from about 10~ 4 if the Higgs mass is just 
above the W + + T threshold to about 10~ 3 if the Higgs mass is just below the tb threshold. 

In a hadron collider, most of the standard production mechanisms for a charged Higgs 
in the mass range below tb threshold involve the production of an additional very massive 
particle ||. The standard production mechanisms for H + include tt production followed by 
the decay t — > H + b, tH + production, W~H + production, and H~H + production. Because 
of the additional very massive particle, events in which a charged Higgs decays into W + T 
will be easily distinguished from W + T events produced by Standard Model processes. 
There is one potentially significant production process for a charged Higgs that could result 
in events that resemble Standard Model W + T events. That process is qb — ► q'bH + , which 
proceeds through a Feynman diagram that involves a virtual W and a virtual top quark || . 

We conclude that it should be possible to observe the production of W ± + T and Z° + T 
from Standard Model proceses at the LHC. The clean experimental signatures for these 
events also makes them valuable lamposts under which to search for heavy particles that 
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decay into W ± + bb or Z° + bb. 
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